FACILITY FORM 60* 


1 


pN> 

|; - a ■ 


'/ - 


IX. 


A 


N 65 -'2 7618 


(ACCESSION NUMBER) 


7^_ 


(NASA CR OR TMX OR AD NUMBER) 


/ 


ICODEJ^ 

Mr 


X-650-b'5-24 


A TA>X- 


DIGITAL COMPUTER SIMULATION 

OF THE 

INFRARED INTERFEROMETER 
SPECTROMETER (IRIS) 

AND 

INTERFEROGRAM ANALYSIS 


GPO PRICE $ . 
OTS PRICE(S) $. 

Hard copy (HC) . 
Microfiche (MF) 




BY 

L. H. BYRNE 


3 


/ 






\ 


JANUARY 1965 


y 


NASA 


GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 



DIGITAL COMPUTER SIMULATION 


OF THE 

infrared interferometer SPECTROMETER (IRIS) 

MD 


IN TERFEROQRAM ANALYSIS 


by 

L. H. Byrne 


GODDARD SPACE FLIGHT CENTER 
GRESNBSLT, MD. 


JANUARY 1 965 



ACKNOWLEDGEMENT 


The author is indebted to the IRIS experimenter. Dr, 

R, Hanel, without whom this work would not have been poss- 
ible. Also, gratitude is due Mr. F. Bartko, Dr. C. Hammer, 
and Dr. R. Roper for the many stimulating discussions on 
the subject topic. Mr, G. Nichols has provided the necess- 
ary information on the electronics of the system, and Mr. 

A. Simmons is credited with the computer coding and his 
tireless assistance in setting up computer and plotter runs. 



A-3-S -T -R-A-C -T 


A digital computer program for the simulation of the 
output of the Infrared Interferometer Spectrometer (IRIS) 
and a program for the analysis of the resultant interfero- 
gram are being effectively utilized for two purposes: (i) de- 
termination and verification of instrument design parameters, 
and (ii) provision of the basis for operational computer pro- 
grams which will be needed to process and analyze the data 
from the actual instrument, 'Both programs are required to 
satisfy (i), and, with proper program design, (ii) is ful- 
filled automatically. The two programs, then, form a "closed- 
loop" synthesis/analysis model for testing the effects of in- 
strumentation and design parameters on both the interferogram 
and the transformed spectral profile. 

Results indicate that a sampling rate for IRIS of approx- 
imately four times the highest frequency present is desirable. 
The effect on the transformed spectral profile due to digitiz- 
ing is, as anticipated, about the same as the effect due to 
the mean expected signal-to-noise ratio at the detector. The 
effect due to finiteness of the optical solid angle is less 
than that caused by digitization or noise at the detector,, 
Additional distortions of the data from the actual interfer- 
ometer may be anticipated due to secondary effects not included 


in the simulation 
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1 .0 INTRODUCTION 


The Infrared Interferometer Spectrometer (IRIS) is a 
Michaels on-type interferometer designed for spacecraft use; 
hence the basic optical principles employed are classical and 
well-known. A beam splitter is exposed to incoming radiation, 
and every radiative component is divided into two separate 
beams. These beams are directed along different paths of the 
optical system before being focused concurrently on the de- 
tector. The optical distance of one of the patns is constant; 
the other path is varied at a constant rate by means of a mov- 
able mirror. The resultant interference pattern, then, is a 
function of path difference, and the measured flux at the de- 
tector is known as the interferogram. 

One of the (advantages of the interferometer is the wide 
spectral range which can be obtained ;\ however a distinct dis- 
advantage is the fact that the spectral analysis must be per- 
formed in terms of the transformed interferogram. This is 
clear when it is recognized that every point of the sampled 
interferogram contains an amplitude contribution due to eacii 
spectral component involved. Thus, the study of many of the 
effects of varying design parameters, as well as the opera- 
tional use of the data, must necessarily be oriented around 
the reconstructed spectral profile. VThe large number of cal- 
culations involved in this transformation process for a wide 
spectral range, and the desired spectral resolution, then. 
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make the use of a fast, large-scale computer a necessity. 

The literature on the theory of interferometric spec- 
tros.copy is abundant, W-W« however most of the work in 
this field was done prior to the advent of the second-gener- 
ation computers of today. As a result, many of the experi- 
menters could not avail themselves of the advantages inherent 
in such a system as described here. 

j The purpose of the programs presented here is two-fold: 

'(i) to provide a means of determining and verifying design 
parameters for IRIS; and (ii) to provide the basis for oper- 
ational computer programs which will be needed to process and 
analyze the data from the actual interferometer. To satisfy 

(i) , a "closed-loop feedback" system, consisting of (a) a sim- 
ulation model of IRIS, and (b) a program to transform the re- 
sulting data, are required. If handled properly, requirement 

(ii) is then fulfilled automatically. 

Any simulation model is, at best, only an approximation 
of certain events. However, these events are usually combined 
in complex fashion, even if the physical phenomenon to be sim- 
ulated appears simple when viewed in terms of effect only. In 
most instances, it suffices to consider only those character- 
istics which most significantly affect the final result. Nev- 
ertheless, use of the model and iteration through the learning 
cycle usually indicate ways in which the model can be improved. 
The "closed-loop" method of synthes is /analysis then also be- 
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comes a means of refining and enhancing the model to whatever 
degree is desirable, as well as serving as an effective tool 
in the initial design stages. 

The most difficult interferometer effect to simulate is 
the total instrument response, since it is made up of a com- 
posite of a number of individual effects. Conceivably, this 
response could even vary from instrument to instrument. The 
individual effects included in the IRIS simulation model are 
those due to (i) the finite solid angle of admission ; (ii) 
the finite path difference (obtained as a direct result of the 
finite mirror travel distance); (iii) noise at the detector; 
and (iv) the digitizing process. The latter, however, is 
external to the instrument itself, and is not included in the 
instrument response. Nevertheless, digitizing has a very de- 
cided effect on the transformed spectral profile 0 

The programs were designed such that they could be used 
for the study of similar interferometer systems merely by con- 
trol of the input parameters; however to avoid over-generali- 
zation, thus preserving simplicity, as well as computer mem- 
ory, certain options and decisions have been left to the dis- 
cretion of the user. The structure of the programs is modular 
and the various effects may be included, or omitted, at the 
user's option, by proper setting of program switches. This 
permits flexible use of the system as a research tool, and min- 
imizes changes required in the event that program modification 
is desiredp 



2.0 MATHEMATICAL ANALYSIS 


The analysis presented herein is intended for use in 
connection with the paper by Hanel and Chaney regarding 
IRIS, Ref, M- The presentation begins with a discussion 
of the "pure" interferogram, i„e., one without any instru- 
mental effects which cause distortion of the data, followed 
by the inclusion of each of the distortion effects, 

2.1 THE 11 PURS" INTSRFEROGRAM AND ITS SiTITHESIS 

The derivation of the equation for the "pure” interfer- 
ogram may -be readily understood by reference to Fig. 1, which 
is a simplified schematic of the interferometer optical system. 



Fig. 1 

Consider an incoming radiation beam of wave number, v Q . 
This beam is divided, by' the beam splitter, into two beams, 


y\j and y^# Bearn is reflected to the movable mirror, , 
and by back to the beam splitter, where it is again divid- 
ed into two components. One of these components is reflected 
by the beam splitter back in the direction of the source, and 
the other is transmitted through the beam splitter toward the 
detector. Beam y^ follows a path to the fixed mirror, M 2 * and 
is also reflected back to the beam splitter, where it, also, 
is divided into two components. One of the latter is trans- 
mitted back toward the source, the other is reflected toward 
the detector. Interference is produced by the two beams arriv- 
ing at the detector due to their respective optical path dis- 
tances; when the traversed distances are equal, the beams are 
in phase. Thus, the state of interference at the detector is 
a fuction of path difference between the two interfering beams, 
and consequently, of the position of the movable mirror, • 

The superposition of two beams of constant amplitudes, 
a.j and a 2 » and the same angular frequency, w^, yields a result- 
ant of constant amplitude, A, and constant phase, ?• [5j- To 
see this, let y^ and y ? be given by 
y,(t) = a sin(w t -<3(0, and 

11°1 (1 ) 

y o (t) = a sin(w t -*(o)» 

where a^,^ and are the amplitude and phase, re- 
spectively, for y^ and Then, the resultant, y, is 

y = y t + y 2 

= a.j sin(w Q t -oC^) + a, sin(w Q t -°C,) 
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(2) 


= (a 1 cosfc.j) + a 2 cos(oC 2 )) sin(w Q t) 

- (a 1 sinfy) + a £ sin$* 2 )) cos(w Q t) 

= A sin(w Q t - <p ) , where 
A cos (f) COS ft, ) -f 3.^ COS fcf?))* and 

t * ~ 

A si n(<p) = (a 1 sinfc^ ) + a £ sin^)). 

A can be found by squaring and adding equations (3) : 
2 2 2 

A = + a 2 + 2a^a 9 cos^ - <y 2 ), 

and <p is found from 


(3) 

w 


ta.n(<p) = (a 1 sinfy) + a 2 sin^ 2 ))/(a 1 cosKj) + a? cos fr 2 ^* ^ 
Thus, for a given w q , it is seen that A and ^ of equation (2) 
are constant, and can be found from equation (A) and (5)* re- 


spectively. 


If the incoming beam is split in half, i.e,, if a^ = a 2 = a, 
equation (A) becomes 

A 2 = 2a 2 (1 + cos£q (6) 

Since the intensity, I , of the beam is proportional to the 


square of its amplitude, vie can write 

I* = A 2 = (B/2)(1 + cos ft, -pf 2 )), (7) 

where B = Aa , 


The difference in phase of the two split beams, however, de- 
pends upon path difference,^", and wave number, v q . Conse- 
quently, equation (7) can be expressed as 

l'(v o ,<0 = (B/2)(1 + cos(2irv 0 «f)), (8) 

where B/2 is the intensity of the original beam. 

Extending (3) over all positive wave numbers, v, and integrat- 
ing with respect to v, equation (8) becomes 
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B(v,T)( 1 + cos(2»veJ)) dv. 


(9) 


iV) 




where B is a function of wave n 'amber, v, and temperature ,T, 

rco 

Hereafter, the constant component, j\ 3(v,T) dv, is disre- 
garded, and equation (9) is taken to be, to : 


I(cf) 


i oo 

B(v, 

0 


T) cos(2irv<() dv. 


( 10 ) 


Equation (10) is the equation of the "pure" interfero- 
gram, where 3(v,T) is the spectral profile, and <f is defined 
over the range (- 00, 00), Taking the range of path differ- 
ence to be infinite is, of course, impossible, since the mov- 
able mirror must traverse a finite distance, (-x^, x^) ; there- 
fore, correspondingly, the effective path difference ,<f, is 
restricted to the range (-Xj, x^). (The effect of this fin- 
ite path difference limitation is discussed in 2,3.2.) Then, 
as<Tvaries from -x^ to x^, time, t, varies from -2^/2 to2f/2, 
and the distance traveled by the movable mirror is 2x^ (note 
that Xj = 2^), Thus, the effective optical path difference 
can be expressed as 

Vjt = 2V m t, (11) 

where V T and V are the constant velocities of the image 
and the mirror, respectively. Also, 

<f ma x = h ?/2 = V m -r - 2x m = Xj, -r/2 i t if/ 2. (12) 

Equation (10), then, becomes 


I(t) 


■r 


B(v,T) cos(4wvV m t) dv, -T/2i:ti ^/2. 


(13) 



How, let 3(v i ,T k ) be the function defined by 
, , ’'S >in 

3(v i>Tk ) = ^ v ] ilin — v ±£ m 

v max^ v i* 

where i = 0,1,2,**», N, and B(v^,T k ) is obtained for 
each temperature, T k , k = 1,2,»»*,K, by means of the Planck 
Blackbody Function: 

B ( v i* T k) = c 6 V i^ exp ip2 V i/ T k^ ” 1 )” 1 » ^5) 

-12 2 -1 

where C, = 2hC„ = 1.1909 X 10 watt cm ster. ; and 
o 2 

Cg = hC/L = 1.^380 cm °K; h is Planck's constant, L the Boltz- 

man constant, C the velocity of light, T the temperature in 

k 

degrees Kelvin, and v. the wave number in cm • 

In order that the bias due to the instrument temperature be 
removed, the spectral profile is taken to be 

3 <VV = ¥W * WV- (161 

where is the blackbody intensity due to the instrument 
temperature, and is the target intensity. . 

Equation (1 3) can be approximated by the sum 

I(tj) = AV 3(V it T k ) 003(4*^7^), (17) 


xjhere v. = v + i(^v), for each 

t . = j(At), j = - M,- M + 1 »***»0,1 
J 

Since (17) is symmetric about t^ = 0, l(t .) may be obtained 
by computing (17) for j = 0,1,2,»»»,M, and reflecting I(t^), 

for all j ^ 0, to obtain l(t .)♦ If, also, the time origin 

"J 

is translated such that the interferogram is defined over the 
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interval t = (0 ,''£), that is, if we let t* = t 4-T/2, or 

t =t' -?/ 2, then equation (17) becomes 

I(t^) a Av B(v i§ T k ) cosC^irv^Ct^ -f/2)), or 

I(t j) = Av ± B(v^,T k ) cos(-4'irv i V^t ) cos(%v^V_'T/2) 


T=0 


+ Av B(v i ,T k ) sinC^^t') sinC^v^^T^) . (18) 


Equation (18) is now the equation of the "pure" digital in- 
terferogram. However, for the purpose of detailed analysis 
and comparison with the transformed spectral profile below, it 
is helpful to carry the mathematical representation one step 
further. 

The maximum resolution obtainable, due to the finite path 
difference, EU [5]. (see, also, 2. 3, 2.2) is 
(Av)' = 1/2XJ = 1/4^, (19) 

and the resolving power, R, obtainable for a given wave num- 
ber, v q , is 

R = v o /(Av)\ (20) 

Thus, using equations (11), (12), and (19)» the velocity of 

the mirror, V , can be written as 
m 

V m = 2x m^ = = l/ 2 ( Av )' • ( 2 0 

Substituting (21) into(l8), the equation of the interferogram 

becomes 

I(t') =AV B(v ,T, ) cos(2nrv j/(Av) '(2M+1)) cosfrv /(A. v) ') 

0 x^O 1 K 1 1 


+ Av 



B(Vi,T k ) 3in(2m i j/(Av) , (2M+1)) sinfrv^/U v)'), (22) 
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Clearly, when v^ /<M) ' is a positive integer, or zero, the sine 
terms vanish and (22) reduces to 

I(t') = 1 Av V B(v. ,T, ) cos(2ffV. j/&v)’(2M+1)). (23) 

w i=o” 1 

This occurs when (a.v) * divides both v min and i(Av) » since 
v i = v min + i & v )« It will become apparent in the sequel that 
this Diophantine relationship is important in the computational 
process, and further, that it imposes a restriction on the use 
of the simulation model and the subsequent analysis of the 
synthesized data. 

2.2 ANALYSIS OF THE IMTERFBROGRAM 

The spectral profile, B(v,T), is the Fourier cosine 
transform of ("10): 

f oo 

I(<f) ccs(2jrv<0 d<f. (24) 

.00 


Taking <i*over the finite range (-x^, x^.) , and substituting 

£= 2V t, and d<f= 2V dt into (24), we obtain 
m m 


B(v,T) = 4V 
m 



l(t) 003(4^7^) dt. 


(25) 


Substituting (21) into (25), applying the time-origin trans- 
formation, t = t' -1T/2, and writing the integrals as summa- 
tions, equation (25) becomes 


B ( v i» T l c ) = (2/((av)’(2M+1))) cos^/Uv)') 



+ (2/((av)'( 2M+1))) sin(nv i / (&v) *) 



I(tj) sin(X), (23) 


where the argument, X = 2»rv i j/((Av) '(2M+1 ))• 
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Equation (26), thsn, is the Fourier transform of equation (22), 
and its form is adaptable to digital algorithmic solution. 

For the simulated data, assuming no phase shift, it 
would be sufficient to analyze only half the interferogram, 
using the cosine transformation. However, in practice, it 
is difficult to determine exact zero path difference, because 
the sampled center point will probably not occur exactly at 
£ = 0. Further, it is possible that the interferogram is not 
symmetric about zero oath difference due to optical misalign- 
ment, and for other reasons. It is, therefore, desirable to 
use a sine and cosine analysis over all the 2M+1 points. 

Then, the spectral amplitude, B C (v^,Tj c ) , may be cal- 
culated from 

B °(via» T k ) = (af a + b i a )% w here (27) 

2 $[ 

a ± = 2 /(&v)*( 2M+1)) 27 I(t!) cos(aH a j/(2M+1)) # and (28) 

a j=0 - 

b i a = 2 /(C*v)*( 2M+1)) £: I(t!) sin(2iri a j/ (2M+1 ) ) , (29) 

j=0 0 

and where i & = v. K&) 

Equations (28) and (29) hold only for integral values of i a # 

For some purposes, it is convenient to observe the 
transformed spectral profile in terms of temperature, T, , as 

K 

a function of wave number, v i a » To accomplish this, it is nec- 
essary to calculate the inverse Planck blackbody values. How- 
ever, 3(vi a ,T k ) may be positive or negative, according as i & 
is even or odd, and as the original 3(v. ,T, ) of the synthesis 





is positive or negative. Since the cosine coefficients pre- 
serve this polarity, it is sufficient to carry the sign of 
a- over to B c . Then, the correctly-signed spectral values 
may be expressed as 

(SIGNUM(a ia )) B c (v ia ,T k ). (30) 

The temperature, T k , for each k = 1,2,«»»,K, may now 
be obtained from 

T k = C 6 v ia /(ln(C 2 v^/((SIGNUM(a ia )) B c (v ia ,T k )) + 1)). (31) 

2.3 DISTORTION EFFECTS 

In section 2.1, the "pure" interferogram was developed. 

It is the purpose of this section to discuss those effects 
which produce the principal distortions in the data. 

2,3.1 THE FINITE SOLID ANGLE 

It- can be shown, [2]. K 1. that the total flux on the 

detector due to a monochromatic source of wave number, v . 

* o 

is given by 

= B ( V 0 » T ) smsin(Tiv/rV2ir)/(irv/ii/2Ti)) cos(2irv 0 rf(l-H/^7r)), 

(32) 

where S is the useful aperture area, and £1 is the ad- 
missable solid angle. 

The (sinX)/X factor and the n/4ir term in the cosine argument 
are perturbations caused by the oblique rays due to the fin- 
ite solid angle. The former represents an amplitude modula- 
tion; the latter takes the form of a slight phase shift in 
path difference and amplitude modulation. 
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For the purpose of the IRIS simulation, S was normalized 
to unity, and since £L is small (= Tr/200 steradians) , il/4rr was 
neglected in comparison to unity. Equation (32) then becomes 
l2( v 0 »^ = |Cv 0 .«> B(v 0 ,T), where (33) 

§(v 0 ,4) = (s in Grv 0 <fi}/21T) / 6r 270 ) cos(2irv 0 <£). ( 3 ^) 

Integrating (33) over all desired, finite values of wave num- 
ber, the equation of the interferogram is now 


i,<«n - 


C' 


max 




vp(v,tO B(v,T) cos(2irv<£) dv, 


min 


(35) 

(36) 


where W.tf) = sin(rrv<fii/27I)/(irv«£iy2'n) , 

The effect of the finite solid angle ,Q, on the transformed 
spectral profile is discussed in sections 2. 3. 2.1 and 2, 3. 2, 2* 


2,3.2 THE LIMITING FUNCTION FOR FINITE PATH DIFFERENCE 

As previously discussed (section 2.1), the restriction 
of mirror travel to the finite distance, (-x^.x^), restricts 
O to the finite range (-Xj, x^). In order to examine the 
effect on the transformed spectral profile due to this finite- 
ness, it is convenient to define the rectangular function, 

D(«f). in such a way that, for a single wave number, v q , equa- 
tion (35) takes the form 

Vv«f) = ^o’f) B (V T > (37) 

where is the rectangular function defined by 

(1 , - x-tf** £ x T 

D(<f) = ) T T (38) 

(_0, elsewhere. 

The effect of this limiting function on the transformed spec- 
tral profile is discussed in sections 2«3»2.1 and 2.3.2. 2. 
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2. 3. 2.1 THE INSTRUMENT FUNCTION 

By definition, the instrument function is the Fourier 
transform of the interferogram resulting from a monochromatic 


ft W. /\^r4- /s in n«+ v»n 1 Vi n vn n* n 

illl CAtCUUCU L,± CLJL X CIU^U ) 


4>'in > /r»v*s noc a cqo no-nr»o 

| >-*. wr w v* V 


of instrument functions of the same form but of varying mag- 


nitudes. 


Consider equation ( 37 ) • Let T (f(x)) denote the Fourier 

c 

cosine transform of f(x), and let f(x) * g(x) be the convo- 


lution of f(x) with g(x). Then, the transformed spectral com- 
ponent, B C (v q ,T), is obtained from 
B C (v o ,T) = T c (I 4 (v o ,*)) 

= T c (DW y(v 0 ,£) B(V 0 ,T) cos(2irv o ^)) 

= T c (D(^)) * T c (f(v 0 ,*)) * T c (B(v 0 ,T) cos(2nv 0 f)) 

= T C (D(*)) * T c (f (v Q ,*)) * T c (l(v 0 ,<f)) 

= T c (Dtf)) * T c («f (v Q( tf)) * T c (T c (B(v 0 ,T))) 

= T C (D(^)) * T c ( iiv 0 ,S)) * B(v q ,T) 

= H(v o ,£) * B(v 0 ,T), ( 39 ) 

where H(v Q ,£) = T q (D(£)) * T c ( ^ (v Q ,$) ) , and l(v,tf) 
is given by equation (10). 

Thus, the transformed spectral component, B C (v Q ,T), is the 

convolution of the observed (source) spectral component, B(v Q ,T), 

and the function, H(v Q ,£). H(v Q ,£) is the instrument function , 

and is the convolution of the (sin X)/X function, T„(D(4)), and 

c 

the rectangular function, T c (y<v« ). It can be seen that 
there is no loss in generality in extending the monochromatic 
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source to the desired finite spectral interval. However, it 
is also apparent that the effect of the internal instrumenta- 
tion on the data must be studied in terms of the instrument 
function for each spectral “component", where the spectral 
"component" is itself a very small spectral interval. 

2. 3,2.2 APODIZATION AMD SPECTRAL RESOLUTION 

The (sin X)/X function, T(D(<T) ) , has a disturbing effect 
on the data for some applications. This is caused by the side 
lobes (secondary maxima) of the (sin X)/X function. Tne pro- 
cess of minimizing this effect is known as apodization . 

It is not the purpose here to give a thorough treatment 
to the topic of apodization— this is a subject for extended 
work on this project. Nevertheless, several significant points 
are in order. 

Apodization is performed at the expense of a loss in 
spectral resolution (2], After linear (viz., triangular) 
apodization, the maximum resolution obtainable, due to the 
finite path difference, is no longer (ref. equation (19)) 

(Av) ' = 1 /2Xj = 1 but rather, 

(Av)' = 1/xj = l/2x m . (40) 

For example, with a mirror travel distance, x^ = 0.2 cm. (as 
is the case for IRIS), (Av) ' = 2.5 cm“^ without apodization, 
and (av) 1 = 5 cm“^ with apodization. 

The finite solid angle ,U, also imposes a limitation on 
the spectral resolution, W. &] . which can be expressed in 
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the form 


( Av ) ' = vJlI/27I 0 (41) 

Again using IRIS as an example, and taking v max = 2000 cm 
andJii = TT/200 ster., it is seen that the maximum resolution 
obtainable for this particular wave number, due to the finite- 

4 

ness of H, is 5 cm” . Therefore, in this case, the limiting 
instrumental resolution factor is the finite solid angle, rath- 
er than the finite path difference. On the other hand, if v 
is taken to be, say, $00 , (&v) * = 5/4 cm”, and the lim- 

iting resolution is no longer due to the finite solid angle, 
2,3*3 N0I5B AT THE DSTECTOR 

The sxgnal-to-noise ratio, S/N, used here, £0 , is the 
ratio of the peak (zero path difference amplitude of the inter- 
ferogram of a blackbody for the maximum expected temperature) 
signal to the RMS value of the noise. The noise is expressed as 


H(Vj) = ^(sJW/RMS), (42) 

where s is a scaling constant; is the peak ampli- 

tude; the x. are normally-distributed (pseudo) random numbers; 

J 


and 


RMS = 


x = 



x) 2 /(2MH 


1 »*. 



x /(2M+1 ) • 
3 


with 


(43) 

(44) 


The equation of the interferogram now takes the form 
I 5 W = I 3 (^) + N(t), (45) 

where I and N are given by equations (35) and (42), 
respectively. 
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2,3.4 DIGITIZING 


;iach interferogram sample of IRIS is digitized into an 

information word of nine bits (eight bits plus a flag bit). 

The dynamic range of the amplitudes in the neighborhood of 

zero path difference is so great, compared to the rest of the 

interferogram, that the equal distribution of all amplitudes 

over this digitizing range would be extremely ineffecient. 

Consequently, a switch in telemetry voltage gain is instituted, 

which, in effect, produces a division of the interferogram 

amplitude by ten when the amplitude exceeds, in absolute value, 

one tenth of the peak value of equation (42) ) , The flag 

bit is used to indicate this switch in gain. 

The net result of this digitizing process is a mapping 

which sends every possible intensity value (+, and 0) onto 

one of 255 + 2(115) = 485 voltage levels, but with a greater 

density of levels (or sensitivity) in the t neighborhood of 

zero intensity. Due to the way in which voltage levels are 

a 

assigned., there are (2 - 1 ) - 485 = 511 - 485 = 26 counts 

which are not utilized out of the possible 51 1 • 

The digitizing system is exemplified by the contents of 
Table 2, 3 . 1 . 1 . 


- 17 - 



3.0 PROGRAM ANALYSIS 


The calculation procedure for each program is discussed 
separately in the following sections: 3«1» The Interferogram 
Synthesis Program; and 3.2, The Interferogram Analysis Pro- 
gram. The computational procedure for each is exemplified by 
the accompanying flow charts, and detailed in the explanations 
of those flow, charts. Program listings are included in Appen- 
dix B. The plot programs are not included here; however plots 
of the results of the sample problem are included in Appendix A. 

The programs are written in FORTRAN II for an IBM-7094 
(Mod. 2)/l401 combination. The core memory capacity of the 
7094 computer used is 32K words. The interferograms are plot- 
ted on an EAI magnetic tape-driven plotter, and the recon- 
structed spectral profiles are plotted on the IBM-1401 printer. 
3.1 THE INTERFEROGRAM SYNTHESIS PROGRAM 

The purpose of the Synthesis Program is to simulate the 
output of the interferometer in the form of both (i) a digi- 
ital sampled interferogram (which, in turn, serves as input 
to the Analysis Program); and (ii) a graphical display (see 6.1) 
of the interferogram. The various effects (ref. section 2.0) 
are applied in the order of their natural occurrence. These 
effects, in their application, are independent of each other, 
and any combination of them, for a given run, is possible by 
means of proper specification of the input parameters. Pro- 
gram notation is defined in Table 1 . 
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TABLE 1 



DSFINITI 

ON 0? NOTATION 

Mathematical 

Notation 

Program 

Notation 

Explanation 

V 

NU 

Wave number 

Av 

DELNU 

Integration interval 

(AV) ' 

DNU 

Resolution 

v min 

NUMIN 

Minimum wave number in spectrum 

V 

max 

NUTOP 

Maximum wave number in spectrum 

v i 

NUONE 

Minimum wave number for each con- 
stant interval of input spectrum 

V 2 

NUTWO 

Maximum wave number for each con- 
stant interval of input spectrum 

i 

I 

Index on wave number 

M 

ITOP 

Maximum value of i 

t* 

TKEP 

Time 

j 

J 

Index on time 

2M+1 

JTOP 

Maximum value of j 

dt' 

DT 

Time increment 

TI 

TI 

Instrument temperature 

TG,T 

TT,TG, TEMP 

Target temperature 

B(v if T) 

BNU(I) 

Blackbody intensity for wave num- 
ber, v^, and temperature, T 

— 

TTNEW 

Target temperature corresponding 
to each interval, (NUONE, NUTWO) 

«*j> 

XINTEN(J) 

Interferogram intensity (amplitude) 
as a function of time samples, t^ 


BTI(l) 

Blackbody intensity due to the 
temperature of the instrument 
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TABLE 1 


DEFINITION OF NOTATION 
(continued) 


Mathematical 

HT L L * 

i UcAb-LUIi 

program 

Notation 

Explanation 

Bra< v i.T) 

BTT(I) 

Blackbody intensity of the 
target 

s 

SCALE 

S/N scaling factor 

N(t‘) 

XNOISE(J) 

Calculated noise as a function 
of time 

V m 

V 

Mirror velocity 

TT 

PI 

3.141 59265 

H 

OMEGA 

Solid angle of admission 

A 

nnax 

PEAK 

Peak signal for noise computation 

* 

XBAR 

Mean value of the noise 

RMS 

RMS 

RMS of the noise 

— 

PEAK2 

Digitizing divide level (=1/10 PEAK) 

AD 

DELD 

Interval between digitizing levels 

— 

SCA(J) 

Digitized (scaled) interferogram 
values 

F 

F 

Transformed spectral interval de- 
sired for Analysis Program 

I 

o 

ZEROIN 

D. C. component (= 0.0) 

a. 

ia 

APRIME(I) 

Transformed cosine coefficients 

b. 

ia 

BPRIME(l) 

Transformed sine coefficients 

K - 1 

KCARDS 

Number of input cards which follow 


the "basic” parameter card 
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1.1 


FL'JA CHART 


CALCULATION PROCEDURE AMD EXPLANATION OF 
The calculation procedure for the Synthesis Program is 
summarized in the form of the following explanation of the 
Flow Chart, 3, 1.1.1, pg. 22, 

BOX 1 

The card input parameters are: SW1 , SW2, SW3» 3N4, NUMIN, 
NUTOP, DELNU , ITOF, JTOP, DT, TI, TT, KCAHDS, and SC. ALE. All 
of these appear on the first card. If KCARDS has a negative, 
or zero, value (see BOX 3 below), this is the only parameter 
cai-d required for one interferogram (a simple, blackbody tar- 
get). If additional spectral structure is required, KCARDS 
must specify the number of spectral components (in addition 
to the basic component specified on the first card) and con- 
sequently, the number of cards to follow the first card. Each 
of these latter cards must contain the target temperature, 
TTNEVf, and the desired wave number interval, (NUONE, NUT'NO), 
for each of the additonal spectral components. (See the sam- 
ple problem for examples of the .card images). 

BOX 2 

Necessary program initialization is performed. 

BOX 3 

The desired spectral profile, B(v^,T^) , is computed in 
accordance with the input parameters and equations (15) and 
(16). If more than a simple blackbody target is desired, one 
of the spectral intervals (usually, the longest, but not nec- 
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START 
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READ 
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Cards 


INITIALIZE 


COMPUTE 

j"Observed" Spectral 
Profile 


COMPUTE 

Right-hand Half of Interf erogram | 
With Finite Solid-Angle Effect 


Yes 
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Set? 


No 
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COMPUTE 

Right-hand Half of 
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REFLECT 

Computed Points To Obtain Left-hand] 
Half of Interf erogram ; Translate 
Time Origin 



0 


APPLY 

Digitizing Effect 
(Scaling) 

9 

Yes / r Sw4 N 

\ Set? / 
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No 
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Program 
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15 
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3. 1.1.1 FLOW CHART. SYNTHESIS PROGRAM 
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essarily) and its corresponding target temperature, must be 
selected as the basic component over the entire spectral range, 
(NUMIN, NUTOP), This, then, is specified aa the first param- 
eter card. The remaining spectral components, determined by 
the contents of the KCAKDS, are superimposed on the basic 
spectrum, prior to the calculation of the composite spectrum, 
by the program. 

BOX k 

SW2 is the "(sin X)/X switch". If S'tJZ is a non-zero, 
positive value, the interferogram is computed with the Finite 
Aperture Effect included (BOX 5b ) . Otherwise, the "pure" in- 
terferogram is computed (BOX 5&)» 

BOX 5a 

The "pure" interferogram is computed in accordance with 
equation (If), for j = 0,1,2,*»»,M. 

BOX 5b 

The right-hand half of the interferogram, including the Fin- 
ite Solid-Angle Effect* is computed (see equation (35)) accord- 
ing to the following: 

I 3 (tj) = Lv £ f<vV B(v i> T k> cos (W/it.), (^) 

where the parameters are as previously defined. 

BOX 6 

The computed points of the right-hand half of the inter- 
ferogram are reflected to produce the entire interferogram, 
which is symmetric about t j = 0 (ref. pg. 8). Then, the time 
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origin is translated to the extreme left, such that t . is de- 

J 

fined over the range (0 X). 

BOX 7 

SW3 is the "noise" switch. If SW3 is set to a non-zero, 
positive value, noise at the detector is calculated and applied 
to each point of the interferogram (BOX 8), Otherwise, the 
program jumps the noise computations and proceeds to BOX 9» 

BOX 8 

Noise is calculated and applied in accordance with the 
procedures noted below and the equations of section 2.3«3» 
Rectangular (pseudo) random numbers, x^, p = 1,2, are 
computed by means of the relationship 

x^ = Kx^ (mod 2 P ), where (47) 

x* is the pth random number; x ,, is the pfl st random 
p 1 p+1 

number; K is a constant (the largest power of five which will 
fit the computer’s word length); and P is the number of bits 
per computer word. 

It is known, £73 » that the set of rectangular random num- 
bers, can be given a near-normal distribution by forming 

the set, ^*x.^, where 


x = 

3 


x*, for each j, 
15=1 P 


(48) 


For this simulation, P* is taken to be 10, Thus, each random 

number, x., of the normally-distributed set, f x -\* i- s sum 
J » J 

of ten rectangular random numbers, x*. The generation of the 

p 

set,£x.|, is accomplished by a FAP subroutine. The remainder 
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of the noise calculations is performed in FORTRAN 


BOX 9 



SW4 is the 

"digitizing switch" 

. If SW4 is a 

non-zero, 

posi 

.tive value, 

the digitizing effe 

ct is applied. 

as described 

below, BOX 10. 

Otherwise, the pros 

ram proceeds directly to 

BOX 

11. 




BOX 

10 





Table 2, below, exemplifies the actual IRIS 

digitizing 

process o Since 

the telemetry data 

format has not 

yet (at the 



TABLE 2 



TM INFO. WORD 




GAIN INFO. 

DECIMAL 

TM 

APPROX * 

BIT 

3ITS 

COUNT 

VOLTAGE 

INTENSITY 

1 

00000000 

0 Not Encoded 


1 

a 

00000001 

1 

+31 .750 

0.?3 X 10- 2 

• 

a 

1 

01110011 

115 

+ 3.200 

0.73 X 10~ 3 

0 

m 

00000001 

1 

+ 3.175 

( ) 

• 

0 

01111111 

12? 

+ 0.025 

( ) 

0 

• 

1 0000000 

128 

0.0 

0.0 

0 

0 

11111110 

255 

- 3.175 

-0.73 x 10 ~ 3 ^ 

1 

• 

10001101 

141 

- 3.200 

( ><h 

* 

0 

1 

11111111 

255 

-31 .750 

-0.73 x 10 “ 2 


DIGITIZING PROCESS 
TABLE 2 


time of this writing) been definitely specified, the mapping 
illustrated by Table 2 is not carried to its final form here. 
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The digitizing effect is applied as detailed below. 

Let I(J) represent the calculated interferogram amplitudes 
(by whatever logical program path this point was reached). Then, 
for all J such that 


I(J)| > j PEAK.2) , for PEAKS = ?EAK,/lO, (4°) 

the I(J) are converted to the new amplitudes, I(J)/10, and those 
points involved are flagged. Thus, PEAKS is the "divide thresh- 
old". A new interferogram, I*(j), results, such that 
I*(J) £ | PEAKS) , for all J. (50) 


It is now necessary to map each of the points, I*(J), onto 
its appropriate digitizing level in the range (-PEAKS, PEAKS). 
This is .accomplished by computing and rounding to the nearest 
integral value the numbers, 3CA(J), where 


SCA(J) = I*(J)/(AD), and where (51) 

AD = PEAKS /1 27. 

Therefore, each of the (scaled) amplitudes, SCA(J), of the 
digitized interferogram is one of the numbers: 

-127, -126,. ..,-1,0, I,***, 126,127. 


BOX 11 

The output tape for printing includes: (i) all input param- 
eters; (ii) the source spectral profile, B(v. ,T. ) ; (lit) noise, 

X K 

if calculated; arnd (iv) the final interferogram, versus time. 

BOX 12 


The resultant interferogram, versus time, is written to 


r 


tape. If the digitizing effect has not been applied, the 
amplitudes are in standard floating point format (the time 
array always is, in any case). If digitizing has been applied, 
the amplitudes are of the form noted in BOX 10. Then, each 
amplitude word is contained in one 7094 36 -bit word according 
to the following format: 

Bit 

Pos. 35 / .../26/25/24/23/22/21/20/19/18/ •••/I/O 
3 ... pxxxxxxxx •••. 

Bit positions 13-25, inclusive, are reserved for the binary 

count designating the digitizing level. Bit position 26 is 

the divide (voltage gain) flag; a 1 in this position indicates 

that the original amplitude was divided by ten. The sign bit 

is carried in bit position 35* Thus, only seven of the eight 

positions, 18 - 25 , are required, but this format will allow 

for easy conversion to the final telemetry word fonnat, when 

required. (Note that the other bit positions are ignored.) 

Intensity and time are written to this tape as separate 
arrays . 

BOX 13a 

The output tape, BOX 12, from the Synthesis Program is 
used as input to the Analysis Program. 

BOX 13b 

The same output tape, BOX 12, from the Synthesis Program 
is used as input to the Interferograra Plot Program. 

BOX 14 

SW1 is used if more than one interferogram is desired in 
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the same computer run. Setting SW1 to a non-zero, positive 
value causes the program to return to BOX 1 and read the next 
set of parameters. Otherwise, the program halts. 


Return to BOX 1 and begin the next interferogram* 

3.2 THIS INTERFEROGRAM ANALYSIS PROGRAM 

The purpose of the Analysis Program is to accept a digital 
interferogram from magnetic tape, transform those amplitudes 
to obtain the calculated spectral profile, B c (v^„,T^), and 
then output the results for printing and plotting. In add- 
ition, since this program and its calculation procedure will 
be used as a basis for the final reduction program, computa- 
tional efficiency is of primary concern (as opposed to the 
Synthesis Program which, because of anticipated limited use, 
employs no special computational algorithms to improve the 
efficiency and running time). The calculation procedure is 
explained in 3*2.1 and exemplified by the Flow Charts, 3.2. 1.1 
and 3.2.1 .2. 

3.2.1 CALCULATION PROCEDURE AMD EXPLANATION OF FLOW CHARTS 
The calculation procedure for the Analysis Program is 
detailed in the following explanation to Flow Charts 3*2.1 .1 
and 3.2. 1.2. The former outlines the gross processing logic, 
and the latter is the detailed computational procedure em- 
ployed in the algorithmic solution which produces the cal- 
ulated spectral profile. Thus, the second Flow Chart is a 
detailed version of part of BOX 5 and al 1 of BOX 6 of 3*2o1«1. 
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3.2.1 .2 FCURI5R THANSF^MAT IOM ALGORITHM FLOW CHART 
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BOX 1 


The input parameters are: 3W1 , SW2, SW3* NUI-UN, NUTOP, 
JTOP, DNU, F , and TJ. (See the sample problem for examples 
of the card Images.) 

BOX 2 

The digital interferogram is read into computer core 
memory from magnetic tape 0 This is the interferogram output 
tape of the Synthesis Program, and contains the time array, 
TKEP(j) , as well as the array, l(j), (or XINTEN(J)). It 
should be noted that the time array, TKEP(J), is not required 
in the Analysis program; however equivalent information is 
required, i.e., the total number of samples, 2M+1 , and the 
time span,^. TKEP(j) is put on the tape for the purpose of 
plotting time versus intensity. Operationally, of course, 
the time of each interferogram must be known in order to com- 
pute satellite position and instrument earth-coverage for 
each interferogram. 

BOX 3 

SVJ2 is the "digitizing descaling switch". That is, if 
the data have had the Digitizing Effect applied (setting of 
SW4 of the Synthesis Program), then SW2 of this program must 
be set to a non-zero, positive value to "descale" the data 
(BOX 4). Otherwise, the program proceeds directly to BOX 5. 
BOX k 

Digitizing descaling consists of reconverting each of 
the numbers, SCA(J), (see 3.1.1, pg. 26) into an intensity 
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value by multiplying SCA(J) by AD, lor each J, and then multi- 
plying each of the flagged (indicating that a switch in gain, 
or divide by ten, was performed on that datum) values by ten. 
This procedure, as well as the conversion back to floating- 
point format, is accomplished by a FAP subroutine, 

BOX 5,6 

BOXES 5 and 6 accomplish the necessary initialisation 
and the sin/cos transformation which results in the calculated 
spectral profile. This procedure is detailed below and exem- 
plified by Flow Chart 3*2.1. 2, 

A solution of equation (23) is sought in the form of 
equations (27), (23), and (29), It can be shown, w. that 
this solution can be computed by means of the following algo- 
rithm: 

a i a = (2 /(Av)'( 2M+1))(I 0 + U u cos(2iri a /(2M+1)) - U 2i ), (52) 

and 

b ±a = (2/(Av)‘(2lvI+1))(U 1i sin(2ni a /(2M+1))), (53) 

where the U.. are defined recursively by 

3 ^* 

U 2M+2,i = U 2M+1,i = 0 ™ 

U ji = I j + 2U j+1fi cos(2iTi a /(2M+1)), for 

j = 2M, 2M-1, •••,1,0, and i = 0,1 ,2, •••.HOP. 

It is now possible to explain the previous statement (2,1, 
pg, 10) concerning the restriction on the use of the simula- 
tion model and the subsequent analysis of the data. Recall 
that (pg, 11) i was defined to be 

3 . 
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i a = ^/(Av)’ = (v min + i (av))/(av)\ f or i = ( 0 , 1 ), and for 
those values of i only that produce integral values of i & » 
The requirement that i be an integer is imposed by the time 
origin transformation, which resulted in equation ( 23 ) » and 
by the Diophantine nature of the algorithm (54 )« Thus, for 
example, if it is desired to analyze the interferogram to a 
resolution of 2„5 cm**' , the integration interval, av, cannot 
be chosen to be, say, Av = 1 cm”'*. For then, i(Av) is not 
•an integral multiple of ^v) * for every i = 0,1, 2, •••,?], In 
this case, the minimum resolution which could be obtained 
from the Analysis Program would be 5 cm“^ • 

It is now convenient to add one more refinement to the 
algorithm. In the event that the user does not wish to re- 


trieve the spectral profile for every spectral component, in 
increments of the minimum resolution, the factor, 'F, is in- 
troduced, F is the " computational resolution", and must sat- 
isfy the condition that it is an integral multiple of (Av)*. 
That is, 

- a = < v ain + 1 (Av)F)/(av) * (55) 

must still be an integer, for all i . The upper limit, ITOP, 
of the range of i (equations (54)), then, is determined by 
ITOP = ( (NUTGP - NUMIM)/F) + 1. (56) 

It should be noted that the program, in iterating on i, 
employs the trigonometric identities for the sine and cosine 
of the sum of two angles, over the range of v = (NUMIN, NUTQP) , 
in increments of F, That is, Av of equation ( 55 ) is taken to 


■ 33 - 


be unity, and the observance of the above restriction is left 
to the discretion of the user. Then, the ith iteration pro- 
duces the sine and cosine of 

2ir(v min + iF)/(Av)'(2M+1), i = 0,1 ,2, • • .,IT0P, 

= (^ain/CAvj'CZH+D) + (2ir£F/(Av) ' (224+1 )) 

= ARG1 + ARG2, 

Therefore, 

cos(ARG1 + ARG2) = cos(ARGI) aos(ARG2) - sin(ARGI) sin(ARG2), 
and 

sin(ARG1 + ARG2) = sin(ARGl) cos(ARG2) + siii(ARG2) cos(ARGl). 

BOX 7 

St/3 is the "amolitu te-to-temperature-conversion switch". 

If 3W3 is set to a non-zero, positive value, the amplitudes 
are converted to blackbody temperatures, T^, Otherwise, the 
program proceeds directly to BOX 9» In either case, the ampli- 
tudes, B c (v^ a ,Ti < ) are included in the printout, BOX 9. 

BOX 8 

The amplitudes, 3IGNUM(a. ) B c (v i ,T, ), are converted to 
blackbody temperatures, T^, by means of the inverse Planck 
Function, equation (31). 

BOX 9 

The output for printing on the IBM-1 ^-01 is written to 
magnetic tape. The printed output consists of: a listing of 
the received interferogram values, I(J), all input parameters 
(BOX 1), I, NU(I), BNU(I) (or B c ), BPRIME(l), APRIME(I), BTT(l) , 
and TMP( I ) (or T). 


- 34 - 


BOX 10 


The arrays, BNU(l) and TSMP(l) , are written to magnetic 
tape to be used as input to the plot program which formats 
the plots of temperature versus wave number for subsequent 
plotting on the 131-1401 printer. 

BOX 11, 12 

SW1 is used if more than one interferogram is desired 
in the same computer run. Setting 5W1 to a non-zero, positive 
value causes the program to return to BOX 1 and read the next 
set of input parameters. Otherwise, the program halts « 

4,0 SAMPLE PROBLEM 

A number of computer nans, on artificial spectra, with 
varying parameters, have been made to date. One of these is 
presented here as an example. Runs were made, using the same 
artificial spectrum, of the following interferograms; (i) "pure 11 
(ii) finite solid-angle effect only; (iii) with noise only; 
(iv) vri-th digitizing effect only; and (v) with all effects. 

All other input parameters were held constant for these five 
cases. The sampling rate used was slightly more than four 
times the highest frequency present (3415 samples for the 10- 
second interferogram). The sixth run was made to include all 
effects at approximately half the sampling rate used above 
(7707 samples). Thus, run (vi) is identical, with respect 
to input parameters, to (v), except for sampling rate. Both 
Synthesis and Analysis runs were made for these six cases. 


1 
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and the plotted output for each is included in Appendix A. 

In addition, an interferogram plot (run (vii)) for a mono- 
chromatic source, v = 65 0, with all effects, is included to 
illustrate the effect of the (sin X)/X faction on the inter- 
ferogram (no analysis output is included for this case). 

The input parameters for run (v) are listed in 6, 1 • 1 
of Appendix A. Except for switch settings, this input is 
also valid for (i) - (iv), inclusive. As noted above, run 
(vi) is identical to (v), except for sampling rate. 

Interferogram plots for the above seven cases are pre- 
sented in sections 6.1.3 - 6.1.9, respectively, of Appendix A. 

The transformed spectral profile, in the form of plots 
of temperature versus wave number , of cases (i) - (vi), in- 
clusive, are shown in 6.1.10 - 6.1.15* respectively. These 
plots are made with a temperature resolution of one degree, 
and a spectral resolution of five wave numbers. Also shown 
on each plot, for comparison purposes, is the original, "ob- 
served" source spectrum. 

In every case, the spectral range is taken over the in- 
terval, v = (251 , 2000) ; the interferogram has a duration of 
10 seconds; the instrument temperature is 270 °K; the input 
spectral profile is unchanged (except for (vii) ) ; and the 
peak signal -to-noise ratio, S/N, is 1423 ; 1 (SCALE = 0.0007). 
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5.0 RESULTS A:!D CONCLUSIONS 


Experience to date with the programs has verified their 
effectiveness in the determination of instrument design and 
in forming the basis for operational data reduction. 

The plot of the "pure" interferogram, 6.1 .3* illustrates 
the large dynamic range in amplitude inherent in the undigi- 
tized interferogram. Concern over this problem led to the 
switch in voltage gain (or "division by ten") and the result- 
ing digitization system previously described herein and in 
The gross effect of the digitizing process on the interfero- 
gram can be seen in 6.1.6; the effect on the transformed spec- 
trum is apparent upon examination of 6.1.13. 

As expected, ; the effect due to the signal-to-noise ratio 
at the detector is about the same as the effect due to digi- 
tizing. The peak signal-to-noise (RMS) initially calculated, 

, was on the order of 1000:1 to 2000:1. Computer runs 
made with varying S/N indicate that 1000:1 is the minimal tol- 
erable threshold. A ratio of 2000:1 would allow excellent re- 
sults in the transformed data. Comparison of the transformed 
spectral profiles for the "pure" case, 6.1.10, and the case 
which includes noise only, 6.1.12, illustrates the effect due 
to noise for this level (as previously noted, 1428:1), 

The effect on the transformed data due to the finite solid 
angle, 6.1.11, appears to be the least of those effects stud- 
ied. However, since the total instrument response in not in- 
cluded in the simulation, more effects on the actual data, in 
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addition to the presently-simulated effect, may be anticipated. 

It may be concluded, then, that any improvement to the 
hardware system, in terms of the above effects, should first 
be made in the signal— to— noise ratio. 

The improvement in the transformed spectral profile due 
to the doubling of the sampling rate has resulted in estab- 
lishing this design criterion at the faster rate, A compari- 
son of the effects of the two rates may be made from the ex- 
amination of 6,1,8 and 6,1,15. In effect, the increase in 
sampling rate results in an extension of the usable spectral 
range by approximately 150-200 cm*’"' , depending upon the S/N, 
and assuming a maximum tolerable error of about i 2-3°K, 
Nevertheless, additional computer runs, utilizing more real- 
istic (and finely-structured) input spectra, should be made 
and analyzed. 

Finally, since the ratio of operational data reduction 
time to real (spacecraft acquisition) time is significant, 
the Analysis Program execution timing is of interest. Of 
course, for either program, execution time is a function of 
a number of variables, hence execution time for only one ex- 
ample is given here. For case (vi), that is, the interfero- 
gram with all effects included for 3^15 samples, the IEM- 709 ^ 
Synthesis Program execution time is approximately 18 minutes. 
For the transformation by the Analysis Program of that same 
interferogram, to a 5 cm resolution across the entire spec- 
tral range, the execution tine is approximately 1,2 minutes 0 
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These times do not include the time for printing or plotting, 
and farther, it must be recognized that additional '•house- 
keeping" programming must be added to the current Analysis 
program in order to produce the operational data reduction 
program. 
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The contents of Appendix B are as follows; 
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INTERFEROGRAM SYNTHESIS. 

SETTING SW1 TO NON-ZERO POSITIVE VALUE CAUSES 

REjAD NEXT INPUT RECORD. 

SETTING SW2 TO NON-ZERO POSITIVE VALUE CAUSES 
INCIDENT RAY EFFECT ( S IN ( X ) /X ) . 

SETTING SW3 TO NON-ZERO POSITIVE VALUE CAUSES 
ADDITION OF NOISE TO SIGNAL. 

SETTING SWA TO NON-ZERO POSITIVE VALUE CAUSES 
DIGITIZING EFFECT. 

STORED PROGRAM CONSTANTS ARE... 

(1) . PLANCK FUNCTION RADIATION CONSTANTS, 

C 6 = 1 . 1909E-1 2 
C2 = 1.4380 

(2) . MIRROR VELOCITY, V = 0.02 CM/SEC. 

(3) . PI = 3.14159265 

(4) . SOLID ANGLE OF ADMISSION, 

OMEGA = PI/200.0 RADIANS. 

(5) . MAXIMUM AMPLITUDE FOR BLACKBODY OF 300 DEGREES K, 

. - PEAK = 0.72 21 0E-02 

D I ME NS I ON TTNEW( 100) ,NUONE( 100 ) ,NUTWO( 100) , 

" " 1TG (2000) » BT I ( 2C00 ) ,BTT(2000 ) ,BNU( 2000) , 

2X1 NTEN( 3425 ) , XNO I S E ( 3425 ) , TKEP ( 342 5 ) , ITI2000) , SC A (3425) 

1 CONTINUE 

READ INPUT TAPE 2 , 2 , S W 1 , SW2 , SW 3 , S W4 , NUM I N , NUTOP , DELNU , 

1 1 TOP, JTOP,DT, TI ,TT ,KCARDS , SCAL E 

2 FORMAT(4F2.0, 214, F3.1, 2 1 4, El 3. 8, 2F3.0, I3.F7.7) 

3 CONTINUE 
XNUMIN .= NUM I N 
V =.02' 

XJTOP = JTOP 
XITOlC= I TOP 

XLMAX = (XJTOP - 1.0) / 2.0 + 1.0 

LM AX = XLMAX 
IF(KCARDS) 7, 7, A 

4 CONTINUE 

DO 6 I = 1 , KC ARDS 

READ INPUT TAPE 2 , 5 , TTNE W ( I ) , N UONE (I) , NUTWO ( I ) 

5 F0RMAT(F3. 0,214) 

6 CONTINUE 
GO TO 9 

7 CONTINUE 

DO 8 I = 1, ITOP 
TGI I ) = TT 

8 CONTINUE 
GO TO 12 

9 DO 10 I = 1 , ITOP 
TG( I ) = TT 

10 CONTINUE 

DO 11 K = 1 , KCARDS 
K1 = NUONE(K) - NUM I N + 1 
K2 = NUTWO ( K ) - NUM I N + 1 
T2 = TTNEW(K) 

DO 11 L = K 1 , K2 
TG ( L ) = T 2 

11 CONTINUE 


6.2.1 SYNTHESIS PROGRAM LISTING 
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INTERFEROGRAM SYNTHESIS. 

12 CONTINUE 

BEGIN PLANCK FUNCTIONS. 

XNU = NUMIN 
C6 = 1.1909E-12 
C2 = 1.4380 

DO 13 I = 1 » I TOP 

EA = EXPFIC2 * XNU/T I ) - 1.0 
BTIII) = C6 * (XNU ** 3.0) /EA 
EB = EXPF ( C2 * XNU/TGl I ) ) - 1.0 
BTT(I) = C6 * (XNU ** 3. 0 ) /EB 
BNU( I ) = BTIII) - BTT ( I ) 

_ ITU! = I + (NUMIN - 1) 

XNU = XNU + DELNU 

13 CONTINUE 

END OF PLANCK FUNCTIONS 

WRITE OUTPUT TAPE 3, 14,NUM IN, NUTOP , DELNU, I TOP, JTOP, DT , TI , TT , 
1SCALE,SW1 ,SW2,SW3,SW4 

14 FORMAT ( 34H1 NUMIN NUTOP DELNU ITOP JTOP , 8X , 2HDT, 10X, 2HTI , 6X , 2 

1HTT,5X,5HSCALE/I6, 5X , I 4, 3X , F4. 1, 2X, 1 4, 2X, 1 4 , 3X ,E 1 3. 8 , 3X , F5. 1 , 3X , F5 
2 . 1 , 2X , F7. 6/ / 26H SW1 S 

3W2 SW3 SW4/lX,F4.1,3X,F4.1,3X,F4.1,3X,F4.1//> 

WRITE OUTPUT TAPu 3 , 140, ( TTNE W( I ) , NUONE ( I ) , NUTWO ( I ) , I =1 , KCARDS ) 
140 FORMAT(20X,8HTTNEW = F4.0, 3X,8HNU0NE = 1 4, 5X , 8HNUTW0 = 14//) 

WRITE OUTPUT TAPE 3,151 

151 FORMAT ( 6H1 NU,10X, 3HBNU,12X, 2HNU, 10X, 3HBNU, 12X, 2HNU, 10X, 3HBNU, 12 

IX, 2HNU, 10 X , 3HBNU , 1 2X , 2HNU , 10X, 3HBNU/ / / ) 

11 = 1 
12 = 5 

_ K9 = IT0P/5 
DO 150 K = 1 , K9 

WRITE OUTPUT TAPE 3, 15 , ( I T ( I ) , BNU( I ) , I = 11,12) 

11=11+5 
12= 12+5 
150 CONTINUE 

15 FORMAT ( 17 ,5X, E12.5,4( 5X, 14, 5X, E12.5) ) 

PI = 3.14159265 
OMEGA = PI/200.0 
F4 = 4.0 * PI * V 
SLOOP = V * OMEGA 
TJ = 0.0 

DO 20 J = L MAX, JTOP 
XNU = XNUMIN 
XINTEN(J) = 0.0 
DO 19 I = 1,1 TOP 
I F ( SW2 ) 17,17,16 

16 SZ = SINF( SLOOP * XNU * TJ)/(SLOOP * XNU * TJ ) 

IF(TJ) 17,17,18 

17 SZ = 1.0 

18 XINTEN(J) = DELNU * SZ • BNU( I )*COSF( F4*XNU • TJ ) + XINTEN(J) 

XNU = XNU + DELNU 

19 CONTINUE 

TJ = TJ + OT 
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INTERFEROGRAM SYNTHESIS 


20 CONTINUE 

BEGIN REFLECTION. 

J2 = JTOP 
J3 = LMAX -1 
00 21 J = 1, J3 
XINTEN(J) = X I NTEN ( J2 ) 

J2 = J2 - 1 

21 CONTINUE 

END REFLECTION. 

PEAK = .722 10E-02 
IF1SW3) 250,250,22 

22 CONTINUE 

CALL RANDOM ( XNO ISE ( 1 ) ,JTOP) 

XE AR = 0.0 

DO 23 J = 1 , J TOP 

XB AR = XB AR + XNOISE(J) 

23 CONTINUE 
SUM = 0.0 

XBAR = XBAR / XJTOP 
DO 24 J = 1 , JTOP 
XNOISE(J) = XNOISE(J) - XBAR 
SUM = SUM + XNOISE(J) ** 2.0 

24 CONTINUE 

RMS = ( SUM/X J TOP ) ** .5 

F NAGLE = ( PE AK * SCALE) / RMS 

DO 25 J = 1 , JTOP 

XNOISE(J) = XNOISE(J) * FNAGLE 

XINTEN(J) = XINTE.N(J) + XNOISE(J) 

25 CONTINUE 
250 CONTINUE 

FORMAT ( 1H1 , 14H RANDOM NOISE.//) 

WRITE OUTPUT TAPE 3,66 

WRITE OUTPUT TAPE 3 , 67 , ( XNO I SE I J ) , J = l.JTOP) 
FORMAT ( IX, 10EL3.5) 

END NOISE COMPUTATION 

WRITE OUTPUT TAPE 3,26 
26' FORMAT (7H1XINTEN/ / ) 

TJ = 0.0 

DO' 27 J = V, JTOP 
TKEP(J) = TJ 
TJ = TJ + DT 

27 CONTINUE 

K4 = JT0P/10 + 1 

K1 = 1 

K2 = 10 

DO 30 K3 = 1 , K4 

WRITE OUTPUT TAPE 3 , 2 8 , ( TKEP ( L ) , L = K1.K2) 
WRITE OUTPUT TAPE 3, 29, ( X I NTEN ( L ) , L = K1.K2) 

28 FORMAT (1X,10E13.5) 

29 FORMAT! 1X.10E13.5//) 
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K1 = K 1 +10 
K2 = K2 + 10 

30 CONTINUE 

IF ( S W 4 ) 36,36,31 

31 CONTINUE 

PEAK2 = PEAK/10.0 
DELO = PEAK2/127.0 

CALL SCAL ( XINTENI 1) , JT0P.PEAK2, DELD.SCAt 1 )) 
WRITE OUTPUT TAPE 3,32 

32 FORMAT! 15H1XINTEN SCALED//) 

K1 = 1 

K2 = 10 

CO 34 K 3 = 1 , K4 

WRITE OUTPUT TAPE 3, 28 , ( TKEP ( L ) , L = K1,K2) 
WRITE OUTPUT TAPE 3 , 29 , ( SCA ( L ) , L = K1.K2) 

33 FORMAT! IX, 10013// ) 

K1 = K1 + 10 

K2 = K2 + 10 

34 CONTINUE 
36 CONTINUE 

CALL WRITEAU INTEN13425) ) 

CALL WRITEA(TKEP( 3425) ) 

IF(SWl) 35,35,1 

35 PAUSE 7 

END(1, 1,0, 0,0, 0,1,1, 0,1, 0,0,0, 0, 0) 
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INTERFEROGRAM ANALYSIS 

SETTING SWI TC NON-ZERO POSITIVE VALUE CAUSES READ NEXT RECORD. 

SETTING SW? TC NON-ZERO POSITIVE VALUE EFFECTS DIGITIZING SCALING. 

SETTING SW 3 TC NON-ZERO POSITIVE VALUE TRANSFORMS A ( I ) • S TO TEMPS. 

DNU IS INTEFEROMETER RESOLUTION. F IS FREQUENCY CONTROL , WHERE 

WAVE NOS. COMPUTED ARE GIVEN BY... 

NU ( I ) = NUMIN + F ( I - 1 ) , FO R I = 1 , 2 , . . . , ( ( NUM AX-NUM I N ) /F ) - 1 . 
STORED PROGRAM CONSTANTS ARE... 

(1) . PI = 3.14159265 

(2) . C-C COMPONENT, ZEROIN = 0.0 

(3) . PLANCK RADIATION FUNCTION CONSTANT, C2 = 1.4380 

(4) . PLANCK RADIATION FUNCTION CONSTANT, C6 = 1.1909E-12 

(5) . CESCALING (DIGITIZING) MULTIPLIER. 

DIMENSION XNU ( 710 ) ,TEMP(710) , X IT ( 7 10 ) , BT I ( 71 0 ) , APR I ME ( 71 0 ) , 

1BPRIME( 710) , IT (7 10 ) , BTT ( 710) ,BNU( 710) , XINTENl 3425 ) 

1 CONTINUE 

READ INPUT TAPE 2 , 2 , SW 1 , SW2 , SW 3, NUMI N , NUTOP , J TOP , DNU , F , T I 

2 FORMAT (3F2.0, 314, F3. 1, F4. 1,F3. 0) 

3 CONTINUE 

CALL RE AD2 ( X INTEN ( 3425 ) ) 

I F ( SW2 ) 7,7,4 

4 CALL DESCAL ( X I NTEN ( 1 ) ) 

WRITE OUTPUT TAPE 3,5 

5 FORMAT ( 1H1, 13H0ESCALED I ( J ) / / ) 

WRITE OUTPUT TAPE 3 , 6 , ( X I NTEN( J ) , J = 1,JT0P) 

6 FORMAT ( IX, 10E13.5 ) 

GO TO 10 

7 CONTINUE 

WRITE OUTPUT TAPE 3,8 

8 FORMAT ( 1H1 , 4H I ( J ) / / ) 

WRITE OUTPUT T AP E 3 , 9 , ( X INT EN ( J ) , J = 1,JT0P) 

9 FORMATt IX, 10E13.5) 

10 CONTINUE 

WRITE OUTPUT TAPE 3,11 

11 FORMAT! 1 HI, 1 1 X , 1HSW1 , 4X, 3HSW2, 4X, 3HSW3.4X, 5HNUM I N , 4X , 5HNUT0P , 4X , 
14H JTOP , 4X , 3HDNU,4X, 1HF, 4X» 2HTI // ) 

WRITE OUTPUT TAPE 3 , 1 2 , S W 1 , SW2 , S W3 , NUM I N , NUTOP , J TOP , DNU , F , T I 

12 FORMAT! F15.0, 2F7.0, 18, 110, I8,F7. 1 ,F6.1,F6.0/ / ) 

C BEGIN FOURIER 

XNUMIN = NUMIN 
XNUTOP = NUTOP 
PM I N = XNUMIN / DNU 

X I TOP = ((XNUTOP - XNUMIN) / F) + 1.0 
ITOP = X I TOP 
X J TOP = JTOP 

PI = 3. 14159265 

ARG1 = 2.0 * PI • PM I N / X JTOP 
ARG2 = 2.0 * PI * F / (DNU * XJTOP) 

C = COSFI ARG1 ) 

S = SINFI ARG1 ) 

Cl = COSF { ARG2 ) 

SI = S 1 NF ( ARG2 ) 

C3 = 2.0 / (DNU * XJTOP) 

DO 14 I = 1 , I TOP 
U2 =0.0 
U1 =0.0 
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J = JTOP 

DO 13 K = 1 , J TOP 

UZERO = XINTENIJ) + 2.0 • U1 * C - U2 
J = J - 1 
U2 = U1 
UI = UZERO 

13 CONTINUE 
ZERO IN = 0.0 

APRIME(I) = C3 » { ZERO IN + Ul • C - U2) 

BPRIMEI I ) = C3 * Ul • S 
Q=C1*C-S1*S 
S=C1*S+S1*C 
C = Q 

14 CONTINUE 

DO 15 I = 1 » I TOP 

BNU(I) = APR I ME ( I ) **2 . 0 + BPRI ME ( I ) **2 . 0 
BNU( I ) = SQRTF ( 8NU ( I ) ) 

15 CONTINUE 

I T { 1 ) = 0 

DO 16 I = 2 * I TOP 

ITU) = 1-1 

16 CONTINUE 

DO 17 I = 1 » I TOP 
X I T I I ) = I T { I ) 

XNUI I ) = XNUMIN + IF * X I T ( I ) ) 

17 CONTINUE 
IFISW3) 24,24,18 

IS CONTINUE 

BEGIN PLANCK 

C6 = l . 1909 E- 1 2 

C2 = 1.4380 

DO 20 I = 1 , I TOP 

E = C2 • XNU ( I ) / T I 

E = EXPF(E) - 1.0 

BT I ( I ) = C6 * I XNU t I ) **3.0) / E 

I F (APR I ME ( I ) ) 19,20,20 

19 BN U ( I ) = -BNU ( I ) 

20 BT T I I ) = BTI { I ) - BNU I I ) 

END PLANCK, BEGIN INVERSE PLANCK 
DO 21 1=1, ITOP 

TEMPI I ) = C6 • XNU ( I ) **3 . 0 / BTT(I) + 1.0 
TEMPI I ) = LOGFITEMPt I ) ) 

TEMPI I ) = C2 * XNUI I ) / TEMPI I ) 

21 CONTINUE 

WRITE OUTPUT TAPE 3,22 

22 FORMAT I 6X , 1HI ,9X, 5HNUI I ) , 1 OX ,6 HBNU I I ) ,7X,9H3PRIME I I ) , 9X , 
19HAPRIMEI I ) ,9X,6HBTT( I ),9X,7HTEMP( I ) ) 

WRITE OUTPUT TAPE 3, 2 3 , { I T I I I , XNU { I ) , BNU I I ) , B PR I ME 1 1 ) , APR I ME I I ) , 
18TTI I ) , TEMP I I ) , I = 1, ITOP) 

23 FORMAT! 1X,I6,F14.1,E18.5,E15.5,E18.5,E16.5,F13.1) 

END OF- INVERSE PLANCK 

CALL WRITE7ITEMPI 710) ) 

CALL WRITE7IXNUI710) ) ^ _y 

24 CONTINUE ' J* ^ / 

WRITE OUTPUT TAPE 3,25 

2 5 FORMAT! 6X , 1H I , 9X , 5HNU I I ) , 1 OX, 6 HBNU I I ) , 7X.9HBPRIME I I ) ,9X, 



INTERFEROGRAM ANALYSIS 


19HAPR IME ( I ) ) 

WRITE OUTPUT TAPE 3 ,26 , ( IT ( I ) ,XNU ( I ) ,BNU( I ) , BPR I ME ( I ) , APR IME ( I ) t 
II = l.ITOP) 

26 FORMAT ( 1X,I6,F14.1,E18.5,E15.5,E18.5) 

IF ( SW 1 ) ^,2^,1 
2g PAUSE 7 
CALL EXIT 

EN D( 1, 1,0, 0,0,0, 1,1,0, 1,0, 0,0, 0,0) 
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00 
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00 
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00 

0 
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00 

0 
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0 
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ENTRY 

WRITE7 

WRITE7 

CLA 

1.4 


STA 

CMD 


WTBA 

7 


RCHA 

CMD ; 


TCO A 

* 


TRC A 

* + 2 


TRA 

CLAM 


BSR A 

7 


TCO A 

• 


TRA 

WRITE 7 

CLAM 

NOP 
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2,4 

CMD 

IORT 

END 
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00000 


ENTRY 

READ2 

00000 

0500 

00 

4 

00001 

READ2 

CLA 

1 » 4 
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0621 

00 

0 

00021 


STA 

CMD 

00002 
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00 

0 
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RTB A 
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00 
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00 
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00 

0 
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00016 

0020 

00 

4 
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2,4 
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2 
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CMD 2 
A 

IORTN 

BSS 
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1 
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3 06541 

0 
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END 
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OOOOO ENTRY 


OOOOO 

0774 

00 

2 

OOOOO 

DESCAL 

AXT 
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0500 

00 

4 

00001 
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00002 

0621 

00 

0 

00036 


$TA 

00003 

0500 

60 

0 

00036 
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-0320 

00 

0 

00037 
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0100 

00 

0 

00026 


TZE 
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60 

0 

00036 

GO 

CAL* 
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-0320 

00 

0 

00043 
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00 

0 

00034 
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0500 

00 

0 

00034 


CLA 

00012 

0771 

00 

0 

00022 


ARS 

00013 

-0501 

00 

0 

00042 
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00 

0 

00042 
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00 

0 
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00 

0 
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60 

0 
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0 
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00 

0 
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60 

0 
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00 

0 

00041 
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60 

0 
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00 

0 
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00 

0 
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NOSCAL 
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00 

0 
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00 

0 
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TRAP 
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1 00001 

2 

0003 2 
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2 
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00 

4 
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00035 
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BSS 
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2 

OOOOO 
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END 


DESCAL 
0,2 
I, A 
CAT 
CAT 

=0000400000000 

NOSCAL 

CAT 

=0400377000000 

TEMP 

TEMP 

18 

=0233000000000 

=0233000000000 

=.56858268E-05 

CAT 

SWIT 

TRAP NON ZERO NO SCALE 

CAT 
= 10.0 
CAT 
TRAP 
SWIT 
GO 

SWIT 
*+ 1 , 2,1 
TERN, 2, 3425 
2,4 
1 
1 

. 0,2 
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